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Abstract 

m 

CN | We review the density matrix formalism and the positivity conditions for general 

multiple spin asymmetries, taking as an example the case p + p— >A + Ain which 
one, two or three spins are analyzed. Some aspects related to quantum information 
and entangled states are discussed. Some positivity domains for pairs and triplets 
of spin parameters are displayed, together with the experimental points. The case 
of inclusive reaction is also treated, taking as an example the spin- and transverse 
Qh! momentum- dependent quark distributions. 
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1. Introduction 



X 

The single- or multiple-spin asymmetries which can be measured using polarized 
beams, polarized targets or analyzing the final-particle polarizations provide important 
information about the elementary processes in particle physics. These spin observables 
may be related by equalities coming from the symmetries of the processes. Besides, they 
satisfy inequalities expressing the positivity of a Grand Density Matrix, R, which describes 
all possible polarized cross sections. The positivity of R insures that the cross section is 
positive for any initial and final spin states, including entangled ones. The resulting in- 
equalities provide consistency checks of experimental data, constrain any parametrization 
of polarized structure functions and, in the future, may be applied to Monte-Carlo event 
generators with spins. 

The positivity conditions are also interesting from the point of view of the quantum 
information carried by spins. The information about the scattering amplitudes is maximal 
when all the independent spin observables are measured. There are inequalities that 
define the allowed domain and are saturated in this case. Conversely, there is a loss of 
information, in other words an increase of entropy, when some particles are not analyzed 
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or in the case of inclusive reactions. In this case, most of the inequalities may be non- 
saturated. 

There is an abundant literature on positivity conditions. Some key papers of the 60 's 
are still relevant. See, for instance, Ref. pQ for a survey and some references. The subject 
has been revisited in recent months due to the results on the reaction 



p + p^A + A, (1.1) 

measured by the PS185 collaboration at CERN [2|. The polarization of the outgoing 
hyperon or antihyperon is revealed by its weak decay. Since a polarized target was used 
in the last runs, observables up to rank 3 can be accessed. 

This contribution is a short introduction to the density matrix formalism and deriva- 
tion of the positivity conditions. In Sections 2-5 we will consider the case of exclusive 
reactions, illutrating it by reaction (1.1). We will also present in Section 6 some re- 
sults of an alternative "empirical" method by which the positivity domains of a subset of 
observables can be very easily discovered. 

Finally, we will consider the case of inclusive reactions and obtain the inequalities 
which must be satisfied by the spin- dependent quark distribution, considered as the prob- 
ability of the elementary splitting process 

nucleon — > quark + X . (1.2) 

for a given longitudinal momentum ratio x = pi(quark)/p^(nucleon). 

2. The spin observables 

The fully polarized differential cross section of (1.1), more generally A + B — > C + D, 
where A, B, C and D have spin ~, can be expressed as 

— (Sa, S B , S c , Sd) = - f — j Cxfiur S£ Sq S t d . (2.1) 

\ / unpol 

The S 's are the polarization vectors of pure spin states (|S| = 1). In the right-hand 
side they are promoted to four- vectors with S° = 1. The indices \,/i,v,t, run from 
to 3, whereas latin indices i, j, k, I, take the values 1, 2, 3, or x, y, z. A summation is 
understood over each repeated index. S x , S y , S z are measured in a triad of unit vectors 
{x, y, z} which may differ from one particle to the other. C\^ VT are the correlation pa- 
rameters. For example, Coooo = 1, C xy oo = A xy is an initial double-spin asymmetry, Coooy 
is the spontaneous polarization of particle D along y, C 0y0y = D yy is a spin transmission 
coefficient from B to D and C^ xy = C xy is a final spin correlation. 

Equation (2.1) also applies to the case of incomplete initial polarizations, replacing 
the unit vector by with |P^| < 1 and the same for B. The final polarizations 
generally depend on the initial ones, e.g., 

Pc s < Sc> = (is)', Vs » as (p "' Pfl - Sc ; Sd = 0) ■ (22) 

\ / unpol 
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3. The density matrix formalism 

Let {|ck)}, with a = ±|, be the basic spin states of particle A, {|/?)} those of B, etc. 
The quantification axis z may differ from one particle to the other. It can be the helicity 
axis p/|p| or the transversity axis, n = x p c / \p A x p c \. We write the spin-dependent 
amplitude of (1.1) as 

( (A, 7 | ® (A, S\ ) M ( \p, a) ® \p, (3) ) = ( 7 5\M\af3) . (3.1) 

For each spin ^ particle we have the single-spin density matrix 

p(P) = i(l + P.<r), (3.2) 

For partial polarizations of A and S, but definite spins of C and .D, the cross-section (2.1) 
becomes 

^ (P A , P b , So, S d ) = trace{ M p(P A ) ® p(P B ) Af+p(S ) ® p(S D )} . (3.3) 

The final two-spin density matrix is 

_ M p{V A )®p{¥ B ) Mt 

" trace{MMt} " (3 ' 4) 

It describes the individual polarizations of C and D and their spin correlations. The 
polarization of the C is obtained by taking the partial trace over the D spin variable : 

pc = tracer p c ,D , i.e., (l\pcW) = ^l(lb\pc,D\l'5) ■ (3.5) 

5 

The lack of information about a system can be measured by various estimators, among 
which the entropy S = — trace{plogp} and the rank of p. Pure states (maximum infor- 
mation) have zero entropy and unit rank. The entropy (resp. rank) of the initial state is 
the sum (resp. product) of the single-particle entropies (resp. ranks). The rank of the final 
density matrix (3.4) is less than or equal to the initial one. Therefore, complete initial 
polarizations ( |P(p)| = |P(p)| = 1) lead to a final pure state. It does not imply that the 
individual polarization of the A, obtained from (3.5), is complete, because the AA state 
may be entangled jSJ. 

Let us now generalize the density matrix in order to describe in an unified way the 
spin correlations inside the final state and the transmission of polarizations (i.e., of spin 
information) between the initial and the final particles. For this purpose we consider the 
fictitious crossed reaction of (1.1), 

(vacuum) — >p + p + A + A . (3.6) 

We restrict this crossing to spin and flavor variables (we do not consider the momenta). 
An initial particle ket becomes a final anti-particle bra of opposite spin, for instance, 

|p,/3> - (p,-(3\ = (p,(3\CPT , (3.7) 
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where C, P and T are the charge conjugation, parity and time- reversal operators (it does 
not matter if reaction (3.6) does not conserve energy- momentum). Accordingly, we can 
rewrite the spin-dependent amplitude (3.1) as 

(7, 8\M\a, (3) = (-a, -(3, 7, 5\ M crossed (vacuum) = (-a, -(3, 7, 5 . (3.8) 

Thus we produce a one-to-one correspondance between the 2-particle transition operator 
M and a 4-particle state vector |\&), which we will call the Grand Wave Function. To 
shorten the equations, we will introduce the notation a = —a, (3 = —(3, etc. For explicit 
values of a, we will use the notations u and d (for "up" and "down", like for quark isospin 
states) instead of +| and — |. Therefore (3.7) will be written 

\u)^(u\, \d)^(d\. (3.9) 

The Grand Density Matrix, R, which describes all possible spin correlations in reaction 
(1.1), is defined by 

(a^8\R\a'^'5') = ( 7 £|M|a/3) (a'fi\M*\i8') = (a^8\V) (p\a'f?i5') . (3.10) 

Like ordinary density matrices, it is hermitian and semi-positive. Its trace is given by 

traced) = (tf = trace(MM f ) , (3.11) 

which is 2 2 times the unpolarized cross section. Dividing by (3.11), we can re-scale R to 
unit trace, as a standard density matrix. As can be seen from (3.10), R describes a pure 
state: R — |^ r )(^ f |, and is therefore of rank one. This is a particular property of exclusive 
reactions. 

The expression (3.4) for the final density matrix can be re-written as 

a trace 5j { R \p l (p) ® p l (p)} } 
P(A ' A) = ^(R) ' (3 ' 12) 

where p t (p) is the transpose of p(p). This transposition is explained in Appendix A. 

The Grand Density Matrix can be expressed in terms of the correlation parameters 
and vice-versa through 

R = 2~ 4 CVr <t\(A) ® a* (B) ® a v {C) ® a T (D) , (3.13) 

Cx^vr = trace{ R [a{(A) (g) (B) ® a v (C) ® a T (D)] } , (3.14) 
where o"o is the unit 2x2 matrix. 

4. Reduction of the density matrix 

It is difficult to have polarized anti-protons. Therefore the practical spin observables 
in reaction (1.1) concern only p, A and A. They are encoded in the sub-density matrix 
i?(p, A, A) = trace s {p*(p) i?(p,p,A, A)} with p(p) = |J, more explicitely 

OMI ^(P> A, A) \pj8') = £<«Wl ^(P> P> A, A) \aP'i5') . (4.1) 
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This density matrix has dimension 8x8, which is still rather large to write down the 
positivity conditions (the original one was 16 x 16). It has a non-zero entropy, brought 
by the p, and rank 2 because a takes two values in (4.1). 

An important simplification occurs due to the symmetry of reaction (1.1) under the 
reflection II about the scattering plane, which reverses the spin components parallel to 
this plane (it is the "B - symmetry" mentionned in [4J. Using the transversity basis, 
where z = n, the scattering matrix M is even under o~ x — > —cr x , o y — > — o~ y , and the 
amplitude (3.1) vanishes when an odd number of transversities are negative. Accordingly, 
the Grand Wave Function \I/ has no components like \uuud) and the original 4-particle 
density matrix R is reduced to 8 x 8. 

For the same reason, the density matrices restricted to fewer particles like -R(p, A, A), 
p(A, A) or p(A) do not mix states with even and odd numbers of ef s. Thus -R(p, A, A) is 
block-diagonal, being the direct sum of two 4x4 matrices, one corresponding to a — u, the 
other to a — din Eq. (4.1). Each of these sub-matrices is of rank one. Similarly, p(A, A) 
is block-diagonal in two 2x2 sub-matrices of rank 2, one corresponding to (a, (3) = (u, u) 
or (d,d), the other to (a, (3) = (u,d) or (d, u) in Eq. (3.12). Finally, p(A) is diagonal, 
which means that the A polarization is normal to the scattering plane. 

5. The positivity conditions 

The positivity of the Grand Density Matrix comes from the very general, but non- 
trivial requirement that the probability of any process is positive. It is not sufficient to 
require that the cross section (2.1) is positive for any set of polarizations {S^, S#, Sc, S^}. 
Let us suppose, for instance, that (2.1) possesses the factor (1 + Sc • So). This factor 
is positive or null for any Sc and S^. However, it corresponds to a final density matrix 
of the form pc,D = [1 + cr c <S> o~ l D }/4: which is non-positive. For example, for the singlet 
spin state, we have o % c ® o~ l D = —3. The probability that reaction (1.1) produces a (A, A) 
system in the singlet state would be negative! Note that the singlet state is entangled. 
This shows that positivity has to be tested with non-entangled and entangled states. 

Similarly, a factor (1 — Sa ■ Sc), which leads to the complete spin reversal Sc = — Sa 
according to (2.2), gives a non-positive R and is therefore forbidden. As an example, let 
us consider the splitting n — ► q + q followed by a quark-hadron scattering q + h — > q' + hi 
where the q is spectator. The intermediate spin correlation is in (1 — S q ■ Sg). If there 
were a complete spin reversal S q = S q > in the quark-hadron scattering, it would lead to 
a final correlation in (1 + S q > ■ Sq), which is forbidden as explained before. 

The general positivity conditions are as follows: a N xN hermitian matrix p is positive 
(respectively semi-positive) if all its eigenvalues are positive (resp. positive or null). 
Let us consider the symmetric functions of the eigenvalues 

Si = r i > S 2 = r i r i ' S 3 = r i r i r k , • ■ • = r x r 2 ■ • ■ r n . (5.2) 

i i<j i<j<k 

S n is the sum of the on-diagonal sub-determinants of order n (when a sub-matrix has 
its diagonal on the diagonal of p, we call it "on-diagonal"). A necessary and sufficient 
condition of positivity, or semi-positivity with N vanishing eigenvalues, is 

S n > for n < N — N , E n = for N — N < n < N . (5.3) 
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If p is (semi-)positive, each of its on-diagonal sub-determinant is (null or) positive. This 
may provide inequalities simpler than, but redundant with (5.3), in the same manner as 
\x 2 \ < 1 is redundant with \x 2 \ + \y 2 \ < 1. 

The matrix p depends on N 2 real parameters. They can be Re(p^/) for i < i' and 
Im(pat) for i < i', or the correlation parameters, which are linear combinations of them. In 
the iV" 2 -dimensional parameter space, the domain of positivity of p is a convex half-cone. 
Its intersection with the hyperplane Ei = trace(p) = 1 is a finite convex domain T>. The 
boundary of V is a sheet of the hypersurface Ejy = det(p) = 0. It is a (N 2 — 2)-dimensional 
manifold of degree N. On this boundary, p is only semi-positive. The other conditions, 
E n > for n = 2, • • • N — 1 define domains which include T>. These "auxiliary" conditions 
serve to eliminate the other sheets of the hypersurface — 0. The hypersurface where 
E n , or any on-diagonal sub-determinant, vanishes is externally tangent to T>. 

As we have seen, for an exclusive reaction the Grand Density Matrix R is of rank 
one. Therefore all the E n 's are vanishing for n > 2 and all the positivity constraints are 
saturated. It can be shown that R is on a "corner" of T>. On the contrary, when much 
information is lost through non-detected or non-analyzed particles, R is "deep inside" T>. 

6. Empirical approach 

The search for inequalities is straightforward using the density matrix method, but 
does not reveal at once the shape of the allowed domains. Also, when one writes the 
conditions on the density matrix, one gets in general a combination of several spin ob- 
servables, and thus one has to make appropriate combinations of inequalities to obtain 
constraints on two or three given observables of interest. 

To circumvent this difficulty, the following method was used in Ref. The real and 
imaginary parts of the amplitudes were chosen randomly, and the spin observables were 
computed using their explicit expression in terms of amplitudes. This detects which pairs 
or triplets of observables fulfill inequalities, and then these inequalities can be derived 
by straightforward calculus. The case of pairs of observables is extensively discussed in 
Ref. and preliminary results on triplets presented at the LEAP2003 conference [§]. A 
sample of the results is displayed in Figs. [T] and EJ 

Notice that only a few types of inequalities are encountered. For pairs of observables, 
say X and Y, each being typically restricted to [—1, +1], one gets the following possibilities 

• nothing: X and Y might reach any point of the square [— 1, +1] 2 , 

• the diskX 2 + F 2 < 1, 

• a triangle AY 2 < (1 + X) 2 . 

For triplets of observables, say X, Y and Z, the following situations are obtained 

• nothing, any point of cube [—1, +1] 3 is allowed, 

• a sphere X 2 + Y 2 + Y 2 < 1, 

• a cone (1 + Z) 2 > AX 2 + AY 2 , 

• a cubic of the type X 2 + Y 2 + Z 2 ± XYZ < 1. 
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Figure 1: Pair of observables restricted to the unit disk (left, here polarisation and 
Cu are shown) or to a triangle (right, here C nn and Cj m are shown). The small grey 
dots correspond to hypothetical, randomly generated, amplitudes, and the larger dots, to 
actual data. 




Figure 2: Triplet of observables restricted to the inner volume a cone or of a cubic. The 
small dots correspond to randomly generated amplitudes, the larger ones (partly hidden) 
to actual data. 

The latter case is the most interesting, since the domain is restricted in space of three 
observables, but each projection cover the whole square, i.e., there is no restriction for 
any pair observables within (X, Y, Z). The border has the shape of a twisted cushion. 
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7. Inclusive case: the spin-dependent parton densities 

As an example of inclusive reaction, let us consider now the elementary process (1.2), 
which we rewrite for fixed momenta and spin vectors as 

N(p,S N )^q(k,S q )+X, (7.1) 

with k = xp + k^. The probability of (7.1) is the spin- and kx- dependent quark density 
in the nucleon. All what we will say below also applies to the quark fragmentation 
q — > baryon + X, only commuting q and N, or to any inclusive reaction of the type 

A] +B ] +X . (7.2) 

already treated by Doncel and Mendez Thus the problem was solved long ago. One 
can also relate the inequalities in (7.2) to those of the crossed reaction [2j 

Ait +A 2 ]^A 3 + X , (7.3) 

by the correspondence S(A 2 ) <-» — S(C). 

For a given spectator state X, the matrix element in spin Hilbert space can be written 
as {(3\Mx\o). Like in (3.6 - 3.7), we consider the fictitious crossed process 

X - A(p, -S N ) + q(k, S q ) , (7.4) 

where \X) = CPT \X). Note that we have also moved the spectator system X to the 
initial state. The Grand Wave Function and Grand Density Matrix are then defined by 

(a(3\y x ) = ((3\M x \a), (7.5) 

fl = £|*x)<*x|. (7-6) 

x 

Thus R corresponds to a statistical mixture. Its rank r is the dimension of the sub-space 
spanned by the vectors \^x) in the (Nq) spin Hilbert space. It cannot exceed the number 
of possible quantum states of the spectator system. In general r > 1, which means that 
some information is lost, taken away by the spectator partons. 

In our case R has dimension 4x4 and depends on a piori 16 correlation parameters 
C^ v through the analog of (3.13). However, like in the 2^2 reaction (1.1), the plane 
defined by p and kr is a symmetry plane. It is therefore convenient to use the transversity 
basis with z normal to this plane, instead of the helicity basis (unless one integrates over 
k<r). In this basis R is even under a x — ► — a x , a y — ► —a y and the only non- vanishing 
coefficients are C 00 = 1, C 0z , C z0 , C zz , C xx , C xy , C yx and C yy . If we take x along p these 
parameters are respectively proportionnal to /i, — f^, hi — hj- T , g±, h± L , gxr and 
h\ + hi T of Ref. |S], all kinematical factors in p^/M or p\j (2M 2 ) removed. However the 
following inequalities are independent of the choice of the x and y axes in the production 
plane. Let us introduce 

1 _ n ± =1= C Z Q _ ± C XX ± Cyy _ ± C X y ± Cy X _ , 

7, = U nn ) o = n n J 7, = U i 7, = V ■ {.<■<) 



s 



Putting the \Nq) basic spin states in the order {\uu), \ud), \du), \dd)}, one has 

(7. 
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+ U + -iV+\ 

o E/- + iv + o 

U~-iV + D~ n + A- 

\u + + iv- o o i>+-4t7 



As expected from the symmetry about the (x, y) plane, R is block-diagonal in two rank-2 
sub-matrices, which obey the separate positivity conditions: 

(Dt) 2 > (^) 2 + (U^ 2 + {V^f , (7.9) 

and 

D± n >0, i.e., \C ZZ \ = \D nn \ < 1, (7.10) 

which agrees with the results of Bacchetta et al jH] and of Ref. jlj. 

If we integrate over k^, the only surviving parameters are Coo = 1, C xx = Aq(x) / q(x) 
and C yy = C zz = 5q(x)/q(x), where q(x), Aq(x) and 5q(x) are the quark number, quark 
helicity and quark transversity 01101 distributions. One obtains the Soffer inequality [TT] : 

25q(x) < q(x) + Aq(x) . (7.11) 

In a simple model of quark distribution where X just consists in a scalar di-quark, 
all the inequalities (7.9 - 7.10) of the k T -dependent case are saturated. Indeed, such an 
object has no spin to carry quantum information away; a fully polarized nucleon delivers 
a fully polarized quark [TP] . This is no more the case if we integrate over the degree of 
freedom k^. Nevertheless, the Soffer bound keeps saturated. 

8. Conclusions and outlook 

The formalism developed in the 60's remains extremely powerful to analyze the con- 
sistency of spin observables. However, it needs some freshening and new methods are 
needed to quickly derive the inequalities within a subset of accessible observables. We 
hope to have worked in this direction. 

One of the basic tool, already used in 0], is a fictitious crossing which gathers all 
the particles on the same side. It is expressed as a partial matrix transposition. Usual 
crossing also links different physical reactions, just reversing the polarization vectors. 

Particle spin physics also touches the more general theory of quantum information, in 
particular with the concept of entanglement. The fact that the particles considered here 
have definite momenta is not essential. The inequalities obtained in particle physics could 
also apply to "gates" between other kinds of quantum information channels like optical 
fibers. 

In a forthcoming article, we shall provide more details about the positivity conditions 
and their physical interpretation. In particular, we will show explicitly that the method 
of the Grand Density Matrix gives the same constraints on observables that the empirical 
approach based on randomly generated amplitudes. 

Appendix A. Effect of crossing on the operators acting on an initial particle 

Together with (3.7), we have (a'\ — > \a'), 

\a)(a'\ -> \a')(a\ , (A.l) 
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and for a linear combination of such elementary operators, 

Y,\oc)A aal (a'\ -> \a')A a a'(a\ , (A.2) 

a, a' a, a' 

Here we have assumed that crossing acts linearly on operators. Indeed (3.7) is the product 
of two anti-linear operations: (i) applying CPT (ii) changing a ket into a bra. Equation 
(A.2) amounts to the matrix transposition A — > A t , provided we choose the same ordering 
for the crossed basis vectors {\u), \d)} as in the initial basis vectors {\u), \d)} (the ordering 
in magnetic number s z is reversed: it becomes {| — |), | + |)})- For a single-spin matrix 
density, the transformation is 

P = ^l + P-v) - f? = pi = ±(l + P'<T), (A.3) 

where P = — P due to spin reversal, and oj_ = —o\ are the Pauli matrices for the {2} 
representation of SU(2) (which is not often used, due to the equivalence {2} <-> {2}). 

Acknowledgements. We thank M. Elchikh for useful correspondence. We also got experi- 
ence in the properties of the positivity domain in collaboration with E. Loyer J2| during 
his academic training at IPNL. 
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